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Abstract  This  report  presents  a  numerical  procedure  for  estimating  the  true  scattering 
coefficient,  from  measurements  made  using  wide-beam  antennas.  The  use  of  wide-beam 
antennas  results  in  an  inaccurate  estimate  of  cr°  if  the  narrow-beam  approximation  is  used  in  the 
retrieval  process  for  <t°.  To  reduce  this  error,  we  propose  a  correction  procedure  that  estimates 
the  error  resulted  from  the  narrow-beam  approximation  and  uses  the  error  to  obtain  a  more 
accurate  estimate  of  a^.  An  exponential  model  has  been  assiuned  to  take  into  account  the 
variation  of  with  incidence  angles,  and  the  model  parameters  are  estimated  from  measured 
data.  Based  on  the  model  and  knowledge  of  the  antenna  pattern,  the  procedure  calculates  the 
error  due  to  the  narrow-beam  i^roximation.  The  procedure  is  shown  to  provide  a  significant 
improvement  in  estimation  of  obtained  with  wide-beam  artennas.  The  proposed  procedure  is 
also  shown  insensitive  to  the  assumed  model . 


1.  Introduction 


Scatterometers,  calibrated  radar  systems,  are  designed  to  estimate,  accurately  and  precisely,  the 
backscattering  coefficient  (r°  of  homogeneous  targets  from  the  received  power.  Because  the 
backscattered  signal  is  the  supeiposition  of  all  scatterers  illuminated  by  the  antenna  beam,  the 
implicit  assumptions  of  this  retrieval  procedure  are  that  those  illuminated  scatterers  have  the 
same  property  as  far  as  backscattering  is  concerned,  and  the  radar  system  has  the  same  response 
to  all  the  backscattered  components.  Practically,  these  asstimptions  are  £q)proximately  valid  only  # 

when  a  narrow-beam  antenna  is  used  in  the  measurement,  and  common  retrieval  procedures  are 
based  on  so-called  narrow-beam  approximation,  which  will  be  explained  in  following  discussion. 

For  a  monostatic  radar,  the  received  power  from  a  distributed  target  is,  mathematically,  related  to 
<r°  by  an  integral  equation,  the  radar  equation,  as  expressed  by  eq.  (1) 


PriO) 


Pr^^Gl  r 
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where 

P,  =  transmitted  power 

X  =  wavelength 

Go  =  maximum  gain  of  antenna 

=  normalized  two-way  anteima  gain  function 
=  scattering  coefficient  of  the  target 
R  -  range  between  antetuia  and  target 
9  -  incidence  angle  of  EM  waves  (see  Fig.  1), 

The  radar  illumination  geometry  is  illustrated  in  Fig.l.  Notice  that  g^,  and  R  are  the 
functions  of  radar  coordinates  {9„,  and  the  integration  is  carried  out  over  the  illuminated  area 
A.  The  received  power  ^  a  function  of  incidence  angle  9,  as  indicated  in  (1). 

Our  objective  is  to  retrieve  from  based  on  (1),  and  assiunptions  are  made  to  facilitate 
the  retrieval  procedure.  One  of  these  assumptions  is  the  narrow-beam  approximation,  which 
states  that,  if  the  antenna  beam  is  sufficiently  narrow,  o’®,  g^,  and  R  can  be  considered 
ipproximately  as  constants  over  the  illuminated  area  A.  Under  the  narrow-beam  approximation, 
variables  cr®,  g^,  and  R  in  (1)  are  replaced  by  their  representative  values  evaluated  at  incidence 
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Fig.  1  lUumination  geometry 
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angle  d,  and  (1)  can  be  simplified  as 


Prie) 


(4;r)3/?<(0) 


(2) 


where  Amid)  =  j  is  the  "effectively"  illuminated  area.  A  common  practice  is  to  define 

A 

Amid)  as  the  area  of  the  effective  antenna  beamwidth  intercepted  by  the  target  plane.  For  a  pencil 
beam,  the  intercepted  area  is  an  ellipse.  Using  (2),  can  be  retrieved  from  measured  data,  radar 
parameters  and  geometry,  # 


a^id)  = 


P,«?)(4/r)V(g) 

P,X^GlAmid) 


(3) 


Equation(3)  gives  so-called  measured  <7°  (denoted  as  ai),  which  is  calculated  under  the  narrow- 
beam  approximation.  Using  (1),  (2)  and  (3),  <T»id)  can  be  related  to  true  as 


<r^id)  *  <T^id) 
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where 

Ridg,^  ) 

Rid„,  d)  =  —  ;■*  **  ,  normalized  range 
Kid) 


*  normalized  scattering  coefficient 


g^ida,  d),  normalized  two-way  antenna  pattern  as  in  (1). 

From  (4)  it  can  be  seen  that  the  o’*  represents  only  an  averaged  value  of  the  true  o®  weighted  by 
the  normalized  antenna  pattern  and  range.  However,  the  condition  of  the  narrow-beam 
approximation  ensures  that  the  integrand  of  (4)  is  approximately  equal  to  one  over  the  integration 
area,  and  <r?,  the  weighted  average  of  o®,  is  close  to  the  true  o’®  at  incidence  angle  d. 

Unfommately,  wide-beam^  antennas  ate  often  used  in  scatterometers  for  various  reasons,  and  the 
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narrow-beam  approximation  is  then  not  a  valid  assumption.  If  eq.  (3)  is  still  used  in  the  retrieval 
process,  <t?  could  deviate  from  the  tme  by  a  large  value.  Ulaby  et  al.  [1,  pp.  755]  briefly 
discusses  the  effect  of  assuming  constant  <r°  in  the  illumination  integral.  With  a  beamwidth  of 
10°,  the  narrow-beam  approximation  could  result  in  an  error  of  several  dBs  in  estimating  (T°  if 
®  true  <T°  changes  rapidly  over  the  illuminated  area. 

Two  examples  of  the  effects  of  the  narrow-beam  approximation  are  shown  in  Fig.2.  The 
theoretical  values  of  of  a  smooth  sea-ice  surface  at  several  angles  are  taken  from  [2,  pp. 
1754],  and  a  smoothly  fitted  curve  is  generated  from  those  data  and  extended  to  the  entire 
•  relevant  range  of  incidence  angles.  The  curve  is  used  to  calculate  received  power  from  (1)  with  a 

pencil-beam  antenna  of  Gaussian  pattern.  The  calculated  power  is  then  used  to  calculate  cr*  from 
(3),  which  is  derived  under  the  narrow-beam  approximation.  As  a  quantitative  illustration,  the 
figure  shows  how  large  the  narrow-beam  qjproximation  error  could  be  if  a  wide-beam  antenna  is 
used  to  measure  0°  when  it  changes  rapidly  over  the  angular  extent  of  the  anterma  beamwidth. 
^  Because  the  of  targets  with  a  smooth  surface  often  exercises  a  sharp  change  for  small 

incidence  angles,  it  is  necessary  to  reduce  the  error  to  an  acceptable  extent  if  the  narrow-beam 
approximation  condition  is  violated  in  the  measurement.  In  other  words,  some  correction  is 
needed  on  o?  to  take  into  account  the  effect  of  average  in  (4),  so  that  the  requisite  accuracy  of 
estimation  is  achieved  even  though  wide-beam  antennas  are  used. 

^  In  section  2,  we  discuss  the  possible  solutions  and  propose  our  approach.  In  sections  3-5,  we 

describe  the  proposed  procedure  in  detail.  In  section  6,  we  present  two  examples,  and  in  section 
7,  we  give  a  summary  along  with  our  conclusion. 


1.  In  this  report,  an  antenna  is  consideied  as  "wide-beam”  if  the  error  in  the  estimate  of  <7°  by  eq.  (3)  is  too  large  to 
be  accepted.  Otherwise,  it  is  narrow-beam.  These  terms  really  involve  three  factors:  antenna  beamwidth, 
variation  of  <7°,  and  error  tolerance  in  estimating  <7”. 
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Fig.  2  Effects  narrow-beam  approximation:  smooth  surface  of  sea  ice 


2.  Correction  of  the  Narrow-beam  Approximation  Error 


*  Given  received  power  {Pr(^*))  measured  at  incidence  angles  {tf*},  we  try  to  estimate  true 

even  though  the  narrow-beam  approximation  may  be  violated.  The  difficulty  of  the 
problem  is  that,  without  the  narrow-beam  approximation,  it  involves  an  integral  equation  (see 
(1)).  Observing  (1),  we  can  consider  the  function  in  the  integrand  as  a 

two-dimensional  spatial  filter.  The  two-dimensional  convolution,  as  represented  by  (1),  states 
^  that  the  input  to  the  filter  is  the  true  cr^,  and  the  output  of  the  filter  is  the  power  received  by  the 

antenna  (notice  that  output  is  one-dimensional,  or  can  be  thought  as  having  rotational  symmetry). 
We  want  to  find  the  input  signal  of  the  filter  from  the  output,  a  deconvolution  problem. 

For  a  general  two-dimensional  deconvolution  problem,  an  inverse  filter  can  be  structured  to 
^  retrieve  the  input  signal  based  on  output  data  as  well  as  characteristics  of  the  filter  and  the 

measurement  noise  [3,  pp.  61].  However,  in  our  case  the  facts  that  the  measurements  are 
generally  made  over  a  small  set  of  incidence  angles  and  that  the  filter  is  space  variant  make  the 
use  of  inverse  filters  questionable. 

Axline  [4]  developed  a  matrix  inversion  procedure  to  invert  the  integral  equation.  One  drawback 
^  of  this  method  is  that  it  needs  Prifik)  to  be  measured  with  small  angular  increments.  If  a 

theoretical  or  empirical  model  can  be  used  to  describe  (7°,  it  is  possible  to  develop  algorithms  to 
estimate  the  parameters  of  the  assumed  model  [5]  [6]. 

Another  approach  for  attacking  the  problem  is,  instead  of  solving  the  integral  equation  (1) 
^  directly,  trying  to  improve  ai,  the  estimate  based  on  the  narrow-beam  q^proximation,  by  adding 

a  correction  term,  in  terms  of  dB,  which  represents  the  error  due  to  the  narrow-beam 
^approximation  in  the  retrieval  process.  In  the  following  discussion,  we  call  the  error  term  a 
narrow-beam  approximation  error,  and  name  this  kind  of  approach  a  narrow-beam 
approximation  plus  correction.  The  narrow-beam  approximation  error  can  be  mathematically 
k  defined  as 


(5) 


\ 


A  mote  accurate  estimate  of  <r®  can  be  obtained  by  adding  this  error  term  to  <r?. 
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where  is  the  estimate  of  E{d^)  and  all  tenns  in  (5)  and  (6)  are  in  dB. 

£(0*)  is  unknown  and  eq.  (5)  cannot  be  used  to  calculate  £(®*)  because  is  unknown. 

However  we  can  assume  that  true  o’®  belongs  to  (or  is  close  to)  a  certain  model  and  the  model 
parameters  can  be  estimated  from  measured  data.  Once  the  model  is  determined,  the  narrow- 
beam  approximation  error  could  be  calculated  via  (4).  The  idea  behind  this  approach  is  that  what 
we  need  is  the  error  ^(0*),  the  difference  between  and  o’®(^*)  in  (4)  due  to  the  effect  of 

wide  antenna  beam,  not  the  model  itself,  so  an  accurate  model  is  not  cracial  as  long  as  the  error 
can  be  estimated  from  the  model  with  accepted  accuracy.  In  this  report,  we  propose  a  numerical 
procedure  to  estimate  the  narrow-beam  approximation  error.  Later  we  will  see  that  the  proposed 
procedure  is  indeed  insensitive  to  assumed  model. 

The  approach  of  narrow-beam  approximation  plus  correction  has  been  used  by  several 
researchers,  e.g.,  Kim  et  al.  [6]  arid  Stiles  et  al.  [7].  Kim  et  al.  proposed  an  exponential  model 
and  a  recursive  algorithm  to  estimate  the  model  parameters.  Stiles  et  al.  described  an  empirical 
formula,  based  on  a  group  of  hypothetical  theoretical  <r®  curves,  to  calculate  the  error  £(0*)  for 
small  incidence  angles.  The  proposed  approach  in  this  report  is  more  general  and  systematic. 


3.  Exponential  Model:  A  Piecewise-Linear  Approximation 


In  order  to  estimate  the  narrow-beam  approximation  error  via  (4),  a  parameterized  <r®  model  as 
the  function  of  incidence  angle  must  be  chosen.  The  choice  of  the  model  is  based  on  the 
following  considerations: 

a.  It  fits  experimental  data,  at  least  approximately. 

b.  It  facilitates  mathematical  handling. 

c.  It  is  robust;  that  is,  deviation  of  true  cr®  from  assumed  model  should  not  cause 
unacceptable  error. 

Note  that  the  purpose  of  introducing  the  model  is  to  estimate  the  error  resulting  from  the  narrow- 
beam  approximation,  so  it  is  possible  to  obtain  a  good  estimate  of  the  error  even  though  the 


model  itself  may  not  accurately  describe  the  true  scattering  characteristics  of  the  target.^ 

A  reasonable  model  is  the  piecewise-exponential  model 

<TV)  =  A,exp(-0/B.)  /  =  (7) 

Generally  N  =  2,  and  the  interval  Ot.)  can  be  identified  by  inspecting  a»  or  "optimum 
fitting." 

There  are  two  reasons  for  choosing  the  exponential  model.  First,  the  exponential  model  has  been 
found  to  fit  reasonably  well  experimental  data  from  North  American  agricultural  terrain,  sea  ice, 
and  ocean[l,  pp.  577].  Second,  the  model  is  simple  and  the  model  parameters  are  easily 
estimated  from  measured  data.  Particularly,  the  two  model  parameters.  A,  and  B,,  are  separable 
so  that  a  combined  approach  of  a  one-dimensional  lookup  table  plus  one-dimensional 
optimization  can  be  established  to  estimate  the  parameters,  as  proposed  in  this  report. 
Mathematically,  the  exponential  model  is  nothing  but  a  piecewise-linear  approximation  of  the 
true  cr°  in  the  log  scale  (dB).  Therefore,  we  can  expect  a  better  model  by  dividing  the  covered 
range  into  more  segments,  and  each  segment  is  modeled  with  distinct  exponential  parameters. 
As  mentioned  earlier,  what  we  are  concerned  with  is  the  narrow-beam  approximation  error  rather 
than  the  model  itself,  and  we  will  give  an  example  to  show  that  even  if  the  trae  a®  deviates  firom 
the  model,  the  results  of  the  procedure  are  still  reasonably  good. 

As  indicated  in  (7),  the  so-called  exponential  model  here  is  actually  a  piecewise-exponential 
model,  or  piecewise-liner  model  in  log  scale;  that  is,  the  model  consists  of  several  segments  and 
each  segment  is  modeled  by  a  distinct  exponential  function.  For  most  measured  data,  at  least  two 
segments  are  needed,  one  covering  small  incidence  angles  and  another  for  large  angles.  Here  the 
model  is  in  contrast  to  the  one  used  by  Kim  et  al.,  who  used  one  exponential  function  to  cover  all 
the  relevant  incidence  angles. 


2.  The  nanow-beam  ^proximadon  enor  is  defined  in  eq.  (5).  and  the  estimated  enor  fiYNn  the  model  is  given  by 
eq.  (18).  Any  modeling  error  will  result  in  errors  in  er%f(.Ot)  (estimate  of  and  in  (estimate  of 

(7?(0k))-  However  if  the  errors  are  positive  correlated,  they  may  catKel  each  other  in  the  estimate  of  the  narrow- 
beam  approximation  error  given  by  eq.  (18) 
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4.  Estimation  of  Model  Parameter  :  Lookup  Table  Method 


4.1  Calculated  cr® 


With  the  chosen  model,  denoted  as  the  so-called  calculated  (denoted  as  <t®)  can  be 
derived  as  (see  (4)) 


r 


8\e,,^a,e) 


<yl{e,^4.)dA 


(8) 


where  Amidk)  is  the  equivalent  area  as  in  (2). 

As  indicated  by  its  name,  cr^  is  obtained  by  numerically  solving  eq.  (8).  On  the  other  hand,  <r»  is 
calculated  from  measured  Pf  based  on  eq.  (3).  The  integration  in  eq.  (4),  which  relates  <t»  to  true 
is  really  embedded  in  the  physical  process  of  the  measurement  although  (4)  and  (8)  have 
exactly  the  same  form.  Another  difference  between  (4)  and  (8)  is  that  tr  0  in  (4)  is  associated 
with  the  measured  target  and  is  unknown,  and  the  parameterized  model  that  we  are  trying 
to  fit  to  measured  data  in  the  hope  that  narrow-beam  approximation  errors  regarding  a*  and  <7° 
are  close  to  each  other.  It  is  very  clear  that  all  other  quantities  appearing  in  (8),  including 
antenna  pattern  and  range,  should  represent  the  same  quantities  involved  in  the  measurement 
process.  Though  errors  are  inevitable  in  modeling  these  quantities,  we  ignore  these  errors  in 
following  discussions.  In  other  words,  values  of  <72  are  experimental  data;  and  on  the  other 
hand,  is  only  a  mathematical  expression  that  reflects  the  effect  of  the  illumination  integral  on  <7° . 

To  avoid  confusion,  some  words  are  needed  here  about  <7°,  a®,  <7?,  and  <7®.  <7®  is  the  true 
backscattering  coefficient,  which  is  never  known  and  which  we  are  trying  to  estimate,  a?  is  the 
measured  backscattering  coefficient  based  on  received  power  and  the  narrow-beam 
approximation  (see  (3));  it  actually  includes  all  effects  of  noise  and  errors  in  the  process  of 
measurement  and  calibration.  <7®  is  the  estimate  of  <7®  based  on  the  assumed  model  and 
estimated  parameters.  The  last  one,  <7®,  the  calculated  backscanering  coefficient,  is  given  by  (8), 
and  it  reflects  the  illumination  integral  process  in  (1),  which  is  ignored  by  the  narrow-beam 
approximation. 
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One  interesting  point  is  that,  even  though  we  are  trying  to  find  the  best  fit  of  measured  data 
within  the  chosen  model,  we  do  not  attempt  to  use  O'®  as  a  final  estimate  of  cr®,  but  only  use  it  to 
find  the  narrow-beam  approximation  error.  This  is  because  we  are  not  sure  whether  or  not  the 
model  (and/or  the  estimated  parameters)  is  accurate.  However  we  believe  that  the  error  estimate 
obtained  using  the  model  is  reasonably  accurate  for  applying  the  correction  procedure.  Fig. 3 
gives  a  concise  illustration  of  the  above  discussion.  Next,  we  discuss  the  parameter  estimation 
procedure  in  detail. 


4.2  Problem  description 


Comparing  (4)  with  (8),  we  can  sec  that  <r?  will  equal  <r®  if  everything  is  a  perfect  match, 
particularly  between  the  parameterized  model  tr®  in  (8)  and  tme  <t®,  and  all  other  error  sources 
are  ignored.  Therefore,  the  difference  between  them  can  be  used  to  assess  the  correctness  of  the 
parameter  estimation.  The  problem  can  be  mathematically  described  as  follows: 

Given  a  measured  vector  <y?  =  (cr?(^|),  •  •  • , where  superscript  t  means  the  transposed 
vector,  over  a  set  of  incidence  angles  [Oj],  we  need  to  find  the  optimum  model  parameters  Aq 
and  Bq  that  minimize  the  squared  Euclidean  distance  between  the  two  vectors,  a,  and  a®, 
defined  by 


=  ll<T?  -  <r®lP  (9) 

=  i(<T.®(^.)-<r®«?,))" 


where  cr®  is  the  so-called  calculated  vector 


whose  components  {(r®(0t))  ^  obtained  from  (8)  with 

<rl{0)=  A  cxpir-e  IB) 

In  (9),  the  error  is  calculated  in  dB  because  <t®  changes  drastically  with  incidence  angles.  In 
addition,  some  kind  of  weighting  can  be  incorporated  into  the  expression  to  en^hasize  pan  of 
the  measured  data. 
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43  Lookup  table  method 


Various  two-dimensional  minimum-finding  procedures  can  be  used  to  find  the  minimum  solution 
of  (9),  but  three  factors  must  be  taken  into  account: 

a.  Does  the  procedure  converge  each  time? 

b.  Does  it  converge  to  the  global  minimum  each  time? 

c.  How  fast  does  it  converge? 

It  is  not  always  easy  to  answer  these  questions.  One  way  to  get  around  the  above  difficulties  is 
the  lookup  table  method.  That  is.  calculate  over  possible  ranges  of  parameters  A  and  B 

and  store  the  results  in  a  table  in  advance.  Then  in  the  retrieval  process,  the  optimum  parameters 
of  the  model  can  be  obtained  by  simply  comparing  measured  data  o*  with  each  entry  of  the  table 
and  finding  the  one  that  has  the  smallest  distance.  Obviously  the  lookup  table  method  can  only 
yield  a  sub-optimum  solution  in  the  sense  that  the  parameters  A  and  B  must  be  quantized  into 
finite  sets  so  that  the  minimum  solution  is  "optimum”  only  as  related  to  the  chosen  parameter 
sets.  The  quantized  error  can  be  reduced  by  using  a  finer  quantization  step,  or  a  larger-size  table. 
On  the  other  hand,  we  have  assumed  that  the  proposed  approach  is  not  model-sensitive  and  we 
are  not  claiming  the  exponential  model  always  accurately  represents  the  true  That  means  the 
narrow-beam  approximation  error  cannot  be  very  sensitive  to  the  model  parameters  either,  so 
that  the  parameters  A  and  B  do  not  need  to  be  quantized  with  very  small  steps.  Otherwise,  the 
whole  procedure  is  questionable.  The  incidence  angle  6  must  also  be  quantized  into  a  finite  set, 
which  is  less  problematic  because  measurements  are  generally  taken  over  a  set  of  pre-designed 
incidence  angles.  Furthermore,  interpolation  techniques  can  be  used  if  ^casionally  the  measured 
data  have  different  angles  from  those  in  the  table.  Another  problem  is  the  memory  requirement 
for  storing  the  table.  Totally  there  are  NgX  N0  values  in  the  table  where  Ng,  N^,  and  Ng 
are  the  sizes  of  quantized  $,  A,  and  B  respectively.  This  requirement  may  impose  serious 
problems  in  practice. 

One  reason  for  us  to  choose  the  exponential  model  is  its  mathematical  advantage.  In  particular, 
the  two  parameters  of  the  model  are  separable,  and  the  "best"  parameter  A  is  easy  to  obtain  if  B 
is  fixed.  As  a  result,  instead  of  a  two-dimensional  table,  only  a  one-dimensional  table  with  B  as 
parameter  is  needed  for  our  model  fitting,  as  seen  soon. 


Observe  that  in  equation  (8),  parameter  A  of  the  exponential  model  only  acts  as  a  scale  factor. 
Therefore,  we  can  separate  the  effects  of  A  and  B  by  expressing  as  A  where  er°(d*) 

is  still  calculated  from  (8)  but  instead  of  using  the  model  in  (7),  the  normalized  model 


=  exp(-0/B) 
is  used  in  the  expression. 

With  the  above  notations,  the  problem  can  be  restated  as 
min  =  min  ll<r?  -  <T?lP 

A,B  A,B 

=  nun  ( n^  ll<r?  -  cr®ll^ ) 

=  min  ( min  lltr?  -  A<r2H^ ) 

B  A 

The  significance  of  (11)  is  that  the  two-dimensional  search  has  been  translated  to  two  one¬ 
dimensional  searches.  This  makes  it  suits^le  to  use  only  a  one-dimensional  lookup  table  to  store 
<t2,  which  is  only  the  function  of  the  model  parameter  B,  and  to  estimate  A  directly  in  the  model¬ 
fitting  process.  Therefore  we  not  only  drastically  reduce  the  storage  requirement  of  the  lookup 
table  method  by  a  factor  the  size  of  quantized  A,  but  also  eliminate  die  quantization  error  of 
A  because  now  it  is  estimated  over  a  continuous  base. 


4.4  Generation  of  lookup  table 


Each  element  of  the  lookup  table  consists  of  a  particular  value  of  B,  say  B,-  out  of  the  set  of 
quantized  |  )  >  ^nd  two  corresponding  vectors,  erf  ,  and  £,  . 

The  vector  erf ,  is  defined  as 

<rf,,  =  (<(^,),--.<yf,(BJ)'  (12) 

where  is  calculated  from  (8)  using  the  normalized  model  in  (10)  with  parameter  Bj. 
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Ei  is  the  error  vector  whose  the  k*’'  element,  is  simply  the  difference  between  and 

exp(-d*  /5,),  the  value  of  the  model  at  0*,  both  in  dB.  This  error  vector  will  be  used  for  the  error 
correction. 

How  to  choose  the  quantized  parameters  set  (B,  ,  i  =  1,2,. ..>1)  depends  on  the  allowable 
memory  size,  the  acceptable  time  spent  on  the  calculation,  and  the  tolerable  error  due  to 
quantization  on  B.  The  quantization  errors  can  be  assessed  by  observing  the  difference  between 
the  components  of  two  vectors^  and  ,  corresponding  to  two  adjacent  quantized  values,  B, 
and  The  set  of  l^.-,  i=l../t)  is  chosen  to  include  all  possible  incidence  angles,  in 
consistence  with  the  designed  experiment. 

For  each  quantized  B,.  is  found  over  the  chosen  incidence  angles  {0j)  by  numerically 

solving  the  integral  (8).  Besides  the  normalized  exponential  model  (10)  with  specified  B,,  the 
normalized  two-way  antenna  pattern  g^(. )  must  be  specified  in  order  to  carry  out  the  integral, 
which  could  be  a  fitted  function  of  measured  antenna  pattern  or  just  measured  antenna  data  if  the 
measurements  were  made  over  fine-enough  angles.  For  the  latter  case,  interpolation  might  be 
needed  in  the  calculation.  The  chosen  integral  area  A  should  be  large  enough  so  diat  the 
contribution  to  the  integral  from  the  outside  of  the  area  can  be  ignored.  Note  that  only  the 
relative  distance  in  (8),  so  the  real  height  of  the  antenna  position  is 

irrelevant  and  any  height  can  be  assumed  for  the  calculation. 

A  Pascal  program  for  generating  a  required  lookup  table  is  listed  in  the  Appendix.  The 
integration  program  used  was  taken  from  [8]  and  slighdy  modified  for  two-dimensional 
integration.  In  the  program,  the  two-way  antenna  gain  was  assumed  to  be  a  Gaussian  pattern. 

Generating  the  lookup  table  is  a  time-consumed  task.  As  a  matter  of  fact,  dw  integration  (8) 
needs  to  be  carried  out  N0  x  Ng  times  where  N,  and  Nf/  are,  respectively,  the  sizes  of  quantized 
0  and  B.  However,  the  table  has  been  made  once  and  for  all  as  long  as  the  same  antenna  is  used. 
The  model  fitting  costs  little  time  once  the  table  has  been  made. 


4J  Parameter  estimation 


The  parameter  estimation  procedure  is  straightforward  once  the  lookup  table  is  made.  Fig.4 
diagrammatically  illustrates  the  procedure.  The  squared  distance  between  a*  and  is  plotted  in 
the  form  of  equi-distance  contours  with  parameters  A  and  B  as  independent  variables.  Parameter 


B  has  been  quantized  into  a  set  of  discrete  values  (5„  i=  and  a  lookup  table  has  been 

made  whose  ith  element  is  a  vector  ,  given  in  (12).  For  each  discrete  B,  (equivalendy,  <t®  ,  , 
the  element  of  the  lookup  table),  the  local  minimum  distance  dJ  is  calcxilated  with  the  fixed  Bj 
and  the  free  parameter  A.  By  comparing  (Df ),  the  model  parameters  are  declared  as  the  ones 
that  yield  the  smallest  distance  D*  among  {Dj).  As  seen  from  the  figure,  quantization  on 
parameter  B  results  in  some  error  in  locating  the  global  minimum.  The  error  can  be  made  as 
small  as  we  wish  by  reducing  the  quantization  step  with  the  price  of  a  larger  table  and  longer 
searching  time.  Also  two-dimensional  inteipolation  techniques  can  be  used  to  reduce  the  error 
further.  The  step-by-step  procedures  for  estimating  the  model  parameters  are  as  follows: 


Step  1  Calculate  vector  from  measured  data  using  (3). 

Step  2  For  each  element  of  the  table,  cr°,  vector,  find  optimum  A,-,  which  minimizes  the 
distance  between  the  measured  vector  <t?  and  =  according  to  the 

following  equation 

A.  =  ^  i(cr^ej)-<rU^j))  (dB)  (13) 

A  ^1 

which  can  be  proven  easily,  as  shown  in  following  analysis. 

For  a  fixed  B,-,  we  want  to  find  the  optimum  parameter  A,-  that  minimizes  the 
squared  distance 

i??  =  ll<T?-A.<T®,ll*  (14) 

=  min  lltr?  -  A<y®  11^ 

A  * 

Notice  that  in  the  above  equation,  tr^s  are  in  dB  as  mentioned  earlier,  Le., 


so 

dJ  =  min  ll<T?(dB)  -  A(d»)  -  <r^(dB)ll^ 


(15) 


Differentiate  (IS)  and  let  the  result  equal  zero,  we  get  (13). 

Step  3  Compare  { )  to  find  the  smallest  one 

D*  mTrdn{Dj}  (16) 

I 

and  the  corresponding  parameters  A*  and  B*,  which  yield  D*,  are  the  fitted 
parameters  of  the  measured  vector  o*  based  on  quantized  parameters  (B,  ).  The 
error  vector  corresponding  to  B*  can  be  read  out  from  the  table  for  error  correction. 


5.  Estimation  of  Backscattering  Coefficient 


Generally,  the  angle  response  of  ai  can  be  divided  into  two  regions  (see  (7)),  by  either  inspecting 
{(7*(B,))  or  trying  different  segmentations  and  finding  a  suitable  division  that  yields  minimum 
combined  D* ,  After  the  model  parameters  for  each  region  ate  determined,  the  corresponding 
error  vectors  can  be  used  to  correct  the  narrow-beam  approximation  errors  on  the  measured  data  # 

( O’*  (^i) ) .  However,  the  correction  at  those  incidence  angles  near  the  intersection  of  the  two  fitted 
model  segments  might  not  be  as  accurate  as  for  other  angles.  A  better  procedure  is  to  use  the 
combined  model  to  directly  calculate  the  error  vector.  With  the  fitted  parameters,  the  combined 
model  is 


A\exp(rOfB]) 
Alexpi-d  I  Bl) 


0^0  <00 
00  ^  0  <  rr/2 


(17) 


The  combined  model  is  again  used  in  (8)  to  calculate  the  error  of  the  narrow-beam 
qrproximation,  £(Bjt) 
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E{.d,)  =  alM-al{0,)  (dB) 


(18) 


where  (T^{0t)  comes  from  (8)  using  the  combined  model  (17).  Because  £{6^)  reflects  the  effea 
of  the  illumination  integral  (1),  we  can  expect  a  better  estimate  of  <t^  to  be 


<t°(tf*)  =  <r2{e,)  +  £(tf*)  idB)  (19) 

We  call  as  corrected  because  unlike  <t?  (uncorrected),  a  correction  factor,  Bid),  has  been 
added  to  overcome  the  error  of  the  narrow-beam  approximation. 


6.  Results 


#  The  proposed  procedure  has  been  tested  on  different  types  of  data,  and  in  this  section  two 
examples  are  given  to  show  die  performance  of  the  procedure.  In  the  first  example,  the  procedure 
was  used  on  data  obtained  by  sampling  the  curves  labeled  "measured  in  Fig.2  from  0**  to  SO** 
with  steps  of  2.5®.  The  corrected  <t®s  ({<^®(^*))),  the  measured  <t®s  ({<r?(^*))),  as  well  as  the  true 
(r®s,  are  plotted  in  Fig.5.  The  figure  clearly  shows  the  improvement  in  the  estimation  of  <7®  using 

•  the  procedure  described.  Notice  that  the  two-segment  exponential  model  is  not  a  suitable  model 
for  the  smooth  surface  <t®  curve,  as  shown  in  the  Fig.2.  The  remarkable  improvement  justifies 
the  earlier  statement  that  the  procedure  is  not  model  sensitive. 

Another  example  uses  true  experimental  data  that  were  collected  by  two  different  systems.^  The 
^  measured  tr®  (cr®)  of  both  systems  and  the  corrected  <7®  (tf®)  for  the  wide-beam  antenna  data  are 

plotted  in  Fig.6.  No  correction  was  made  on  the  narrow-beam  antenna  data  because,  except  at 


3.  Two  radar  systems  were  used  to  record  data  simultaneously  during  the  CRREL  1989  e;q)eiiment.  One  radar  was 
a  FM-CW  system  with  an  antenna  of  3‘’  beamwidth,  and  another  was  a  step-frequency  system  with  an  antenna 
of  IS**  beamwidth. 
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Fig.  5  Performance  of  correction:  smooth  surface  of  sea  ice 


Fig.  6  Comparison  of  narrow  and  wide  beam  data  (C  band,  1-17-89) 


normal  incidence,  the  narrow-beam  approximation  errors  were  very  small  and  negligible.  The 
figure  shows  that,  as  expected,  the  wide-beam  antenna  system  only  gives  an  erroneous  estimate 
of  cr®  for  small  incidence  angles  if  the  narrow -beam  approximation  is  used  to  retrieve  a®.  On  the 
other  hand,  data  from  both  systems  show  good  consistency  after  the  correction  was  made  on  the 
wide-beam  data. 


7.  Summary 


In  this  report,  we  have  presented  a  correction  procedure  for  the  estimation  of  <t°  for  data 
collected  by  a  wide-beam  system.  The  procedure  estimates  the  error  incurred  by  the  narrow- 
beam  approximation  and  provides  a  more  accurate  estimate  of  <t®  by  correcting  the  error  in  some 
degree.  An  exponential  model  is  used  in  the  procedure,  and  a  combined  ^proach  of  a  one¬ 
dimensional  lookup  table  plus  one-dimensional  optimization  is  proposed  for  model  fitting.  The 
procedure  is  shown  to  provide  significant  improvement  on  the  estimate  of  (r°  and  to  be 
insensitive  to  the  model  assumption  of  <r®. 


» 
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9.  Appendix  Pascal  Program  For  Generating  Lookup  Table 


PRCXJRAM  BEAMINT: 

CONST 

PI  =  3.1415926536: 

HHGHT  =  1; 

H2*l;  |H2-H*H) 

D_R  =  0.0174532925; 

BEAMWIDTH_3DB  =  0.2617994 ;  {  15  DEGREES  5  GHZ  } 

1  BEAMWIDTH_3DB  =  0.0994838  ;  5.7  DEGREES  9.6  GHZ  } 

1  BEAMWIDTH_3DB  *  0.0877028  ;  5.025  DEGREES  13.6  GHZ  ) 

S_THETA  =  2.5;  {  DEGREE  | 

TYPE 

GLNARRAY  »  ARRAY[1..10]  OF  DOUBLE; 

VAR 

THETA  JHUVLFBEAM_INT.CENTER_Y4-ONG_AXIS,SHORT_AXIS, XI,  Y1,Y2,Y0 :  DOUBLE; 
RHOl jmO2J>HIl.SINT.Sl_CIRCLE,S2_CIRClEJl.R0,X_Y :  DOT’^'  . 
R_THETAEAR,S_PARJ14.SIGMA41ALFBEAM_EFF,AILL  :  DOUBLE; 
cmcLEji,;; :  integer: 

FDAT :  TEXT; 


FUNCTION  DB(y:DOUBLE):DOUBLE; 
BEGIN 

DB  :=  lO^LNPO/LNlIO) 

END; 

FUNCTION  TAN(X:DOUBLE):DOUBLE; 
VAR 

C:  DOUBLE; 

BEGIN 
C COS(X); 

IF(C  «  0)  THEN  TAN  1E36 
ELSE  TAN  :-SIN(X)/C 

END; 


{  arccos  function  ) 

FUNCTION  ACOS(C;DOUBLE)-.DOUBLE; 
VAR 

S>VC :  DOUBLE; 


BEGIN 

IF  (ABS(C)>1)  THEN  BEGIN 
IF(C>0)THENC:=1 
ELSE  C:=-l 

END; 

IF(C=0)THEN 
AC  :=  PI/2 
ELSE 
BEGIN 

S:=SQRT(1.0-SQR(C)); 

AC  :=  ARCTAN(S/C); 

IF  (AC  <  0)  THEN  AC :»  PI+AC 
END; 

ACOS :»  AC 
END; 


{  antenna  gain  function  ) 

FUNCTION  ANTENNA_PArrERN(PHI :  DOUBLE):DOUBLE; 
CONST 
0=1.6651092; 

BEGIN 

ANTENNA.PATTERN  :»EXP(-SQR(C*PHI/BEAMWIDTH_3DB)) 
END; 


{  sigmaO  function :  exponential  model  | 

FUNCTION  SIGMAO(PHI.PAR  :  DOUBLE):DOUBLE ; 
BEGIN 

SIGMAO  ;=  EXP(-PHI/PAR) 

END; 


{  find  angle  between  the  antenna  loolt_direction  and  ground  a  pixel  ) 

FUNCTION  ANGLE(RHOFHIJl,YOaiO :  DOUBLE):DOUBLE; 
VAR 

CPm  -.DOUBLE; 

BEGIN 

CPm (RHO*COS(PHI)*YO+H2)/(R*RO); 

ANGLE  :=  ACOS(CPHI) 

END; 


{  ANTENNA  GAIN  fimction  :g**2  } 

FUNCTION  FUNC(PHI,RHO:DOUBLE):DOUBLE; 
VAR 
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ANG  iDOUBLE; 

BEGIN 

ANG :»  ANGLE(RHO J»HI4i.Y0JR0): 
FUNC  :aANTENNA_PArTERN(ANG); 
END; 


{  find  parameters  of  illuminatioa  ellipse  ) 

PROCEDURE  ELLIPSE(THETA,HALFBEAM_INT :  DOUBLE ; 

VAR  CENTER_YJLONG_AXIS,SHORT_AXIS  :  DOUBLE); 

VAR 

A3.C.D£ ;  DOUBLE; 

BEGIN 

A :»  SQR(TAN(HALFBEAM_IND); 

B  :aSQR(COS(THETA))-A*SQR(SINCraErA)); 

C ;»  HEIGHT*(1+A)*SIN(2*THETA); 

D  :=  SQR(HEIGHD*(A*SQR(COS(THETA))-SQR(SIN(THETA))): 
CENTER_Y  :=  C/(2*B); 

E  :=  D+SQR(C)/(4*B); 

LONG.AXIS  :=SQRT(E/B): 

SHORT.AXIS  :=  SQRT(E) 

END; 


(  find  angle  of  a  ellipse  given  range  in  polar  coordinate  system  } 

FUNCTION  PHI_LIMIT(CENTER_Y,L0NG_AXIS,SH0RT_AXIS3H0  :  DOUBLE):DOUBLE; 
CONST 

NOT.ZERO  =  lE-30; 

VAR 

A3,C.D.Y_C0F^_C0F3IMrr :  DOUBLE; 

BEGIN 

Y_COF :«  SQR(LONG_AXIS); 

X_COF SQR(SHORT_AXIS); 

A :»  SQR(RHO)^(Y_COF-X_COF); 

B  :=  2*X_COF*CENTER_Y*RHO; 

C  :=  X_COF*Y_COF-Y_COF*SQR(RHO>X_COF*SQR(CENTER_Y); 

D ;» (.B+SQRT(SQR(BM*A*C))/(2*A); 

LIMIT  :=ACOS(D); 

IF  (LIMIT  =  PHI  1 )  THEN  LIMIT PHIl  +  NOT.ZERO; 

PHI.LIMIT  :=  LIMIT; 

END; 


(polynomial  extrapolation  function  ) 

PROCEDURE  POLINT(XA,YA  :  GLNARRAY;  N  ;  INTEGER  ;X  :  DOUBLE ;  VAR  YJDY  IX)UBLE); 
VAR 
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NS;^4  ;  INTEGER; 
W.HP,HOT)IFT,DIFT)EN ;  DOUBLE; 
CX> :  GLNARRAY; 


BEGIN 

NS:»1; 

DIF:-ABS(X-XA[1]); 

FORl:=lTONDO 

BEGIN 

DIFT:»ABS(X-XA[I]); 
IF  (DIET  <  DIF)  THEN 
BEGIN 
NS:=I; 

DIF:=DIFT; 

END; 

C[I]:*YA[I]: 

D[I];=YAn]; 

END; 

Y:*YA(NS]; 

NS:»NS-1; 

FOR  M:«l  TO  N-1  DO 
BEGIN 

FORI;-lTON-MDO 
BEGIN 
HO:»XAa]*X; 
HP;-XA[1+M1-X; 
W:-Cn+l]-D[I]; 
DEN:=HO-HP; 
IF(DEN-0)THEN 
II :»  I; 

DEN:aW/DEN; 

D[I]:=HP*DEN; 

C[l]:»HO*DEN; 

END; 

IF(2*NS<N-M)THEN 

DY:»C(NS+1I 

ELSE 

BEGIN 

DY:-D[NS1; 

NS:-NS-1; 

END; 

Y:»Y+DY; 

END 

END; 


{ inlegral  along  phi  direction  with  fixed  range  in  polar  coordinate  system  ) 


PROCEDURE  TRAPZD  J»HI(YAB  ;  DOUBLE;  VAR  S :  DOUBLE; 
N  :  INTEGER:  VAR  GUT :  INTEGER); 


VAR 
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J :  INTEGER: 

X,TNM,SUMJDEL :  DOUBLE; 

BEGIN 

IF(N  =  DTHEN 
BEGIN 

S:=0.5*(B-A)*(FUNC(A,Y)+FUNC(B.Y)): 

GUT:=1: 

END 

ELSE 

BEGIN 

TNM:»GLIT, 

DEL:*(B-A)/rNM; 

X:»A40.5*DEL; 

SUM:»0.0; 

FOR  J:=l  TO  GLIT  DO 
BEGIN 

SUM:-SUM+FUNC(X,Y); 

X:-X+DEL: 

END; 

S:*0.5*(S4<B-A)*SUMAmi); 

GUT:»2*GLrr; 

END 

END; 


I  integral  along  phi  direction  with  fixed  lange  in  polar  coordinate  system  ) 

FUNCTION  QR0MB_PHI(YAB  :  DOUBLE):DOUBLE; 

LABEL  99: 

CONST 
EPS  =  l.OE-6: 

JMAX»  16: 

JMAXP-  17: 

K-4: 


VAR 

U.GUT_PHI ;  INTEGER; 

DSS,SS :  DOUBLE; 

H,S  :  ARRAY[L.JMAXP]  OF  DOUBLE; 
CJD  :  GLNARRAY; 

BEGIN 

H[l]:-1.0; 

FORJ:»lTOJMAXDO 

BEGIN 

TRAPZD_PHI(YA3.S[J]J,GLIT_PHI); 

IF(J>«K)THEN 

BEGIN 

FORI;-  ITOKDO 
BEGIN 


END; 

POLINT(CJ).K,0.0,SSJ)SS); 

IF  (ABS(DSS)  <  (EPS*ABS(SS)))  THEN  GOTO  99; 
END; 

S[J+1]:-S[J1; 

H[J+1];»0.25*H[J]; 

END; 

n:*2; 

99:  QROMB_PHI  :*  SS; 

END; 


{ integral  along  range  direction  on  a  function  whose  points  are  phijdirection 
integrals  in  polar  coordinate  system  } 

PROCEDURE  TRAPZD_RH0(A3  :  DOUBLE;  VAR  S  ;  DOUBLE; 

N  ;  INTEGER;  VAR  GUT :  INTEGER;  CIRCLE  :  INTEGER); 

VAR 

J:  INTEGER; 

RHO,TNM,SUMJDELJPHI2_A.PHI2_B,PHI2_RHO :  DOUBLE; 

BEGIN 

IF(N«l)THEN 

BEGIN 

IF  (CIRCLE  o  0 )  THEN  (  integral  on  a  8emi_cirde  | 

BEGIN  (  C3RC1£  >  2  fn  left  semLciicle  ) 

PHI2_A  :=  CIRCLE*PI/2;  (  C3RCLP,  =  1  for  rigltt  semi_circle  ) 
PHI2_B  PHI2_A  (  integral  limits  am  PI/2  to  PI  or  0  to  PI/2  ) 

END 
ELSE 
BEGIN 

PHI2_A  :=*  PHI_LIMIT(CENTER_YJLONG_AXIS,SHORT_AXIS  A); 
PH12_B :»  PHI_lJMIT(CENTER_Y4-ONG_AXIS3HOKr_AXlS3) 
END; 

R :«  SQRT(SQR(A)+SQR(HmGHT)); 

R4:-SQR(SQR(R)); 

SIGMA SIGMA0(ARCTAN(A/HEIGHT)3AR); 
S:«A*QR0MB_PHI(A3HIU»HI2_A)*SIGMA/R4; 

R :»  SQRT(SQR(B>+SQR(HEIGHT)); 

R4  :»SQR(SQR(R)); 

SIGMA SIGMA(KARCTAN(B/HEIGHT)3AR); 

S  :-S+B*QR0MB_PHI(B3HI13HI2_B)*SIGMA/R4; 

S:-0.5*(B-A)*S; 

G1JT:«1; 

END 

ELSE 

BEGIN 

TNMt-GUT; 

DEL:«(B-A)/rNM; 


RHO;»A+0.5*DEL; 

SUM:»0.0; 

FORJ:»lTOGUTDO 

BEGIN 

IF(CIRCLEoO)THEN 

BEGIN 

PHI2_RH0  :=CIRCLE*PI/2: 

END 

ELSE 

BEGIN  (  integral  limit  on  ellipse  edge  } 

PHI2_RHO  :=  PHI_LIMIT(CENTER_YJ,ONG_AXIS,SHORT_AXIS  JIHO); 
END; 

R  ;=  SQRT(SQR(RHO>fSQR(HEIGHD);  (  range  } 

R4 :»  SQR(SQR(R)): 

SIGMA  SIGMA(XARCrAN(RHO/HEIGHT)  J>AR);  {  sigmaO  | 

SUM:=SUM+RHO*QROMB_PHI(RHOJ‘HIlJ»HI2_RHO)*SIGMA/R4: 

RHO:*RHO+DEL; 

END; 

S:=0.5^(S-HB-A)*SUM/rNM); 

GLIT;-2*GLIT; 

END 

END; 


PROCEDURE  QR0MB_RH0(A3  :  DOUBLE;  VAR  SS  :  DOUBLE;  CIRCLE  JNTBGER); 
LABEL  99; 

CONST 
EPS  -  l.OE-S; 

JMAX-16; 

JMAXP-17; 

K»4; 


VAR 

U,GLIT_RHO :  INTEGER; 

DSS  :  DOUBLE; 

H.S  :  ARRAY[1.JMAX?]  OF  DOUBLE; 

CJ) :  GLNARRAY; 

BEGIN 

H[l):-I.O; 

FORJ:»lTOJMAXDO 

BEGIN 

TRAPZD_RH0(A3.SIJU.0Lrr_RH0,CIRCLE); 

IF(J>-K)THEN 

BEGIN 

PORI;«  ITOKDO 
BEGIN 

C[I]  :-H[J-K+Il; 


POLINT(CJ).K,O.O.SSJ>SS); 

IF  (ABS(DSS)  <  (EPS*ABS(SS)))  THEN  GOTO  99; 
END; 

S[J+1];=S[J]; 

H[J+I]:=0.25*H[J]; 

END; 
n  :=*  3; 

99;; 

END; 


BEGIN 

ASSIGN(FDAT.’G5_S.DAT’); 

HALFBEAM.INT  :*  2*BEAMWIDTH_3DB;  ( integnd  beam  ) 

HALPBEAM_EFF  :=  1.2*BEAMWIDTH_3DB/2;  {  effective  beam  ;  illuminadoa  area  ) 
S_PAR  :■  0.05*D_R;  |  exponential  parameter  step  } 

PAR  :s  0.8*D_R;  ( initial  value  of  exponential  para.  } 

WHILE  (PAR  <a  (2.5*D_R))  DO 
BEGIN 

APPEND(FDAr); 

WRITELN(FDAT); 

WRITELN(FDAT4»AR/D_R); 

CLOSECFDAT); 

THETA  ;a  0;  ( initial  inddenoe  angle  ) 

WHILE  (THETA  <51  )DO 
BEGIN 

n  :b  0;  (  n  >■  0  indicate  finishing  integral  correctly  ) 

R_THETA  ;■  THETA*D_R; 

RO  :■  HEIGHT/COS(R_THETA); 

YO :»  HE1GHT*TAN(R_TOETA); 

ELLIPSE(R_THETAJIM.FBEAM_INT.CENTER_Ya.ONG_AXIS;SHORT_AXIS); 
S1_(3RCLE  :»0; 

S2_CIRCLE 0; 

SINT  :»0; 

CIRCLE :»  0; 

PHIl  :=0; 

IF  (CENTER.Y  <  LONG.AXIS)  THEN 
BEGIN 

{ integral  on  left  semi.drcle  | 

CIRCLE :»  2; 

PHIl  PI/2;  ( low  integral  limit  for  angle  ) 

XI  :«  SHORT_AXIS*SQRT(l.SQR(CENTER_Y/LONG_AXIS)): 

Y1  :=L0NG_AXIS-C:ENTER_Y; 

IF(X1  >  Y1)THENX.Y  :aXl  {  up  integral  limit  for  range  ) 

ELSEX_Y:-Yl; 

QROMB_RHO(0.0,X_Y.S2_aR<XB.CIRCLE); 

( integral  on  right  semi.drde  ) 

CIRCLE  ;■  1; 

PHIl  :>  0;  ( low  ituegral  limit  for  angle  ) 

QROMB_RHO(0.0^1^1_CIRCLE.CIRCLE); 
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CIRCLE  :sO; 

RHOl  :=  XI;  {  low  integral  limit  for  range  ) 

RH02  CENTER_Y+LONG_AXIS  {  up  integral  limit  for  range  } 

END 

ELSE 

BEGIN 

RHOl  :*  CENTER_Y-L0NG_AXIS:  |  low  integral  limit  for  range  } 

RH02  ;=  CENTErIy+LONG_AXIS  {  up  integral  limit  for  range  } 

END; 

IF  (RH02  >  RHOl)  THEN 
( integral  on  ellipse  ) 

QR0MB_RH0(RH01.RH02.SINT, CIRCLE); 

SINT  :a  SINT+S1_CIRCLE+S2_CIRCLE; 

ELLIPSE(R_THETAHALFBEAM_EFF,CENTER_YJ.0NG_AXIS,SH0RT_AXIS); 

I  calculate  illuminated  area  ) 

AILL  :=  PI*LONG_AXIS*SHORT_AXIS: 

SINT 2*SINT*SQR(SQR(R0))/AILL; 

WRITELN(’PAR,THErA.SIGMAO^  » ’J»AR:10:5.THErA;10:5,SINT:10;7, 
SINT/SIGMAO(R_THETAJ»AR):  10:7); 

IF(noO)THEN 

WR1TELN( ’ERROR  IN  ’^AR/D_R:10:5,’  WITH  ’JI); 

APPEND(FDAT); 

WRrrELN(FDArj)B(SINT):  10;6JDB(SINT/SIGMA0(R_THETAJ»AR)):  10:6); 
CLOSE(FDAT); 

FOR  JJ  :=  1  TO  5  DO 
WRITELN; 

THETA  :=  THETA  +  S_THETA 
END; 

PAR  :»PAR+S_PAR 
END 


